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\ Abstract 

, In this Letter we propose a multi-Higgs extension of the standard model with 

' Abelian and non-Abehan discrete symmetries in which the mass matrices of the 

charged fermions obtained from renormahzable interactions are diagonal. However, 
(— I I non-diagonal contributions, that are important for obtaining the CKM matrix in 

Qh' the quark sector, arise from non-renormalizable dimension five interactions. Active 

p) • neutrinos acquire mass only from non-renormalizable interactions, the non-diagonal 

O , entries arising through dimension five operators, while the diagonal entries comes 

from dimension six operators. Realistic mixing matrices in the neutrino and the 
quarks sectors are obtained. 
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The flavor problem can be divided in several subproblems |1]: Why do weak 
isospin partners have different masses? Why are quark and lepton masses split? 
^ I Why is there a mass hierarchy between generations? Presently we know that 

' neutrino oscillation data are well described by massive neutrinos p] and the 

flavor problem becomes more interesting: why is there a mixing angle hierar- 
chy in the quark sector but not in the lepton sector? Usually mass matrices 
have the form Ma/3 = J2i(Xi)ai3{^i), where the Fj's are, for Dirac fermions, 
arbitrary complex dimensionless 3x3 matrices, and ($°) denotes one or a 
set of vacuum expectation values (VEVs) of the neutral scalar fields in the 
model. For Majorana fermions, the Fj's are complex symmetric matrices. The 
mixing matrix and the mass pattern of each charged sector depend on the 
structure of the respective Fj's. It is well known that explicit and predictive 
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forms of these matrices can be obtained by imposing flavor symmetries. For 
instance, it has been considered global unitary (horizontal) symmetries like 
[/(l)'s, SU{2)h [3] or SU{3)h More recently symmetries like Sn, [5j; 
A4 0, ^'^^T [IIH], the double covering of A4; Qq, the double dihedral group [Sj, 
and other discrete symmetries have been also considered. 

Here we will turn the problem upside down. Mass matrices will be of the form 
Ma/3 = /(('^*°))a/3 in which / deuotcs a few dimensionless parameters 
(one at leading order), and ($°) is a matrix built with the VEVs of several 
scalar fields. At first sight, there is no gain in predictive power, we are just 
changing a dimensionless general matrix Ta/B by another one with mass di- 
mension (($°))o/3. However, it seems easier, at least in principle, to explain 
patterns of dynamical variables like VEVs, than dimensionless numbers. The 
value of the former can be explained by the dynamics (for instance by study- 
ing the scalar potential) and extra flavor symmetries that we can impose to 
the model. 

In this Letter we will consider an extension of the standard model in which 
the charged fermion masses, that arise from renormalizable interactions, are 
diagonal. However, dimension 5 operators [TT| give important non-diagonal 
contributions to quarks but not to the charged leptons. On the other hand, the 
neutrinos mass matrix is generated mainly from dimension 5 and 6 operators, 
with the former having zero trace as in the Zee's mode [12]. The dimension 6 
operators induce non-negligible diagonal terms and a general 3x3 symmetric 
mass matrix for those neutrinos is obtained. The non-diagonal elements are 
important in the quark sector because they involve also the VEV of the singlet 
which it is not constrained by the electroweak scale, and also in the neutrino 
sector because theses entries involve ratios of doublets' VEVs that produce an 
enhancement. This enhancement does not occurs in the charged lepton sector 
because in this case the non-diagonal entries depend on the same doublet's 
VEVs that appear in the diagonal ones but are suppressed by A~^. No right- 
handed neutrinos are introduced. In summary, we show that all the questions 
mentioned above have a natural answer: fermion masses and mixing arise from 
different scalar sectors and from different effective interactions. The predictive 
power is a consequence of the discrete symmetries imposed to the model. 

Simple non-Abelian discrete groups have been considered for the first time 
as candidates for flavor symmetries by Frampton and Kephart [7j. Later, ^4 
symmetry has been used to obtain realistic mixing matrices in the lepton 
sector [13], and then also in the quark sectors [2]. More recently [U], by using 
these non-Abelian discrete symmetries it was possible to obtain a tribimaximal 
neutrino mixing matrix [16]. The basic idea is to consider SU{2)l fermion 
doublets and singlets in the triplet representation of the ^4 symmetry and 
Higgs multiplets transforming as triplet or singlet of A4. 
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Let us briefly review some properties of the representations. Denoting the 
vector basis (3) as (oi, 02, a^) and (bi, 62, ^3) we have 3®3 = 3^©3b©1© 
1'© 1". Exphcitly dH, 



[3 (g) 3]a = Sa = (a2&3, ash, ai&2), [3 © 3]b = 3b = (a3&2, aibs, a2&i), 
[3 © 3]i = aibi + a2&2 + 0363, [3 © 3]i' = aibi + uja2b2 + uj'^asbs, 
[3 © 3]i" = ai6i + co'^a2&2 + t^ct3&3, (1) 

and 1 © 1 = 1, 1' © 1" = 1, 1' © 1' = 1", 1" © 1" = 1'. 

We consider a model with G^m © © Z3 © Z3 © Z'^ symmetry, where Gsm is 
the gauge symmetry SU{3)c © SU{2)i © f/(l)y. The matter flelds transform 
as usual under Gsm, left-handed lepton and quark doublets L^, a = e,fi,T; 
Qi, i = 1,2, 3, and their respective right-handed singlets laR, um and dm. In 
the scalar sector, we introduce scalar doublets with Y = +1: Hi, Hi, H[, iff, 
and x; three scalar triplets with Y = +2, denoted by % [18], and 
flnally one complex (with Y = 0) scalar singlet (. Fermion weak eigenstates 
flelds transform under {A4, Z3, Z'^, Z'^) as follows: 



L = (Le, Lf,, Lr) ~ (3, u, io^, 1); Ir = {fXR, tr, cr) ~ (3, 1, 1, u), 

Ql = (QiL, Q2L, Qsl) ~ (3, uj"^, 1, 1); ur= {cr, tR, ur) ~ (3, 1, 1, u); 

dR = {sR,bR,dR) ^ {3,1, uj,uj), (2) 

in which co = e^'^*/^. Notice that all fermion flelds are still symmetry eigen- 
states but we do not use an special notation. On the other hand, scalar flelds 
transform under (^4,^3,^3, Z'^) as: 



H = {H,,H2,H,) ~ {3,u,uj^,l), H = {H^,H2,H^) ~ {3,uj,uj\uj^), 
H' = {H[,H'^,H'^) ~ (3,u,l,uj), H" = {H'I,HIH'^ ~ {3,u'' 
$ = ($1, $2, $3) ~ (3, 1, 1, 1), $ ' = ($'1, %) ~ (3, u, CO, io), 

T = (Ti, T2, Ts) ~ (3, u;^ ^, 1), C ~ (1, 1, ^, !)• (3) 



With the flelds in ([2]) and ([2]) we have the leading contributions to the lepton 
and quark Yukawa interactions: 



^=j(^^ [{L-eH)]A[{Le^)]B + ^[L'ea ■ f]B[L^%x + ..•)+ H.c. 
+ (gi[LH]A + 4^[LH]a\C\' + --]ir + H.c. 
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+ (^KIQlH'U + ^[Ql<^']bC + --^un + B.C. 
+ (h,[QlH"]A + ^[QI<^"]bC + --]dR + H.c, 



(4) 



where A is an energy scale characterizing an unknown physics; and e = ia2 
((72 being the Pauh matrix), [ ]a,b means the appropriate product defined in 
Eq. ([1]) and ■ ■ ■ denotes higher- dimensional operators, we assume that \ V\ <^ A, 
where V denotes the VEV of any of the scalar of the model. 

The contributions to the quark mass matrices obtained from renormalizable 
interactions are diagonal, but the non-renormalizable dimension 5 interactions 
contribute with important and not negligible non-diagonal elements, on the 
other hand, for the charged leptons this corrections come from dimension 6 
operators are negligible. We obtain diagonal entries to the mass matrices as a 
consequence of the definition of the basis in which the product [[ ]a-]i gives the 
diagonal contribution while the product [[ gives the non-diagonal contri- 
butions to the mass matrix. For instance, for the u-quarks, the leading diago- 
nal mass matrix arises through the interaction [QH']a = {Q2H^, Q^H[, QiH!^) 
and then [[QH']A{cR,tR,UR)]i = Q2H'r^CR + Q^H[tR + QiH'^UR. Similarly, using 
the [ ]b product, we obtain that [[<5$']b(cr, t/?, ti/?)]! = Q^^^^cr + Qi^'^tR + 
Q2^'iUr, which gives the non-diagonal contributions to the charged lepton 
mass matrix. 



Hence, from the Yukawa interactions in Eq. 
are 



the mass matrices obtained 



9i 



^ V2 ^ ^ 



{)3 
{;i 



B.C.. 



(5) 



for the charged leptons. We have denoted (ffj) = Vi. Notice that the renor- 
malizable interactions are the dominant one, i.e., g[V\ ~ m^-, giVj, ~ and 
giV2 — vHf,- In the quark sector we have 



V'2 






















+ B.C., 



(6) 
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and 







+ H.c, 



(7) 



in which = and ad = jf"-^, for the 2/3 and —1/3 charged quarks, 

respectively. We have denoted {Hf) = Vi, {H'^) = v[, {H'^^) = v-, (<l>i) = v^^, 
($^) = v'^., = f^., and (C) = V(^. Notice, as well, that each charged sector 
has its private VEVs, thus all these matrices are independent from each other. 
In the case of neutrinos we have the mass matrix 






V2 








V2 





VS ^<t>2 


Vl 


Vl v^^ 


1^2. 


VS '"'t'2 






^ 6i ^ 



^2 
(53 



A 



where the 6i are given by (5i = ^ {v^^v^^v^/ Aviv^J, ^2 = iv^^v^/ Kvi), and 
= jt K2^02^xM^i^0i)' where (x) = and (A°) = v^^. 

It is interesting that the discrete symmetries of the model implies that the 
charged lepton mass matrix, Eq. ([5]), is practically diagonal, and that the 
quark mass matrices have three zeros, see Eqs. ^ and ([7]). To the best of 
our knowledge, this mass matrices have not been considered in literature |19j . 
Notice that in the case of the neutrino we have already taken into account the 
hermitian conjugate. 

Since the only constraint on the VEVs of the doublets is that "^^ < (174GeV)^ 
(there are triplets as well, see below), a possibility is that only v'-^ < (174 GeV)^, 
while all the other VEVs are <C 174 GeV. In this case v[ give the main contri- 
bution to the mass of the t quark. On the other hand, the sum of the VEVs 
of the triplets has to have a upper limit of the order of a few GeV to not spoil 
the observed value of the p parameter [20] . Hence, |^|/A are well-defined ex- 
pansion parameters. The present model is natural in the Yukawa interactions 
i.e., unlike the standard model, there is no a hierarchy among the Yukawa 
couplings. In practice, in the analysis below we will consider that all Yukawa 
coupling are in the interval 0.1 — 3. The VEV of the singlet is only con- 
strained by the condition vc^/K < 1. Below we will use always, just as an 
illustration: A = 1 TeV and vc_ = 140 GeV. However, recall that the singlet ( 
may be related to a new energy scale different from the electroweak scale. 

As a first illustration, let us consider the matrix in Eq. ([H]). We will as- 
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sume, without losing generality, that this matrix is a traceless matrix i.e., 
that the condition = ^ is valid [I2] . This is obtained by imposing (using 
Vx = VxG^^'' 1 where and 9^ are real): 



^T3 + ^<t>Z 


— 01 — 6^ 


>1 


TT, e^^ - 




= 0; 9r^2 + ^(Pi 


(1/2) 


fu 






fu 








u 










fu Aviv^^ 





Hence, the neutrino mass matrix is of the form pT][2^ 



^ -26 sine cosO^ 



M" ^ mo 



sin 6 6 

cos 9 e 



e 
6 



(9) 



(10) 



where mo = {fu/A)\v^^vi\, \v2/vi \ = sm9, \v^Jv^^\ = cos 9, e = \{v3V^^/viv^^)\. 
From Table I of Ref. [22] we can see that neutrinos masses and mixing an- 
gles are in agreement with all experimental data if, for instance, in Eq. (fTOj) . 
ml = Aml^ = Amlt^ = (0.039 eV)^, e = 0.21, 6 = 0.236, and 9 = 0.84. In our 
model, these values for the above parameters are obtained by assuming: A = 1 
TeV, f 1 ~ 1 MeV, f<^j = Vi^^, 2vt^ = 2vt2 = "^Ts- with vt^ = 1-25 GeV; then 
we obtain: f^v^, ~ 40 keV, v^^^ = AO/f^ keV, V2 = 0.1 keV, V3 = 0.2 MeV, 
ifu/fl)vx = 0.2 GeV. We also obtain with these VEVs, Amj^ = 9 x 10-^eV^ 
sin^ 29q = 0.82, sin^ 29atm = 0.99, and sin 6^13 = 0.019 a value that is in agree- 
ment with the CHOOZ limit Ves = sin ^13 < 0.07 [23j. We see that a realistic 
leptonic mixing matrix arises since the mass matrix of the charged leptons 
is almost diagonal. If fu, fl ~ C'(l), it implies that some VEVs are of the 
order of MeV or even less depending on the values of and A. This does 
not necessarily implies the existence of light scalars since for having them in a 
model, there must be some symmetries that allow them to appear in the mass 
spectra: The symmetries in the present model do not play that role and /i^^- 
terms survive in the scalar mass matrices (see Ref. [21]). 

On the other hand, with M"-'^ given by Eq. §i) and ([7]) we see that M"('^)M"('^)t ^ 
]Vf"(<^)t]Vf"('^). Hence we need four matrices, V^'^ , for diagonalizing these mass 

matrices by biunitary transformations and the constraint Vckm = ^l^l'^ ■ 
Besides, the matrices V^'^ must satisfy the constraints VfM'^M'^^V^^^ = 
diag{m%ml,ml) = {M^f and Vj^ W^Vj^^ = {M'^f; similarly for the 
V^^ matrices. The quark masses and the Vckm matrix arise, for instance, us- 
ing for the d-type quarks = 0.14, h^i = h'^ = 0.1, and (all VEVs and masses 
below are in GeV) v'l = 42, v!{ = 1.4, w^' = 0.062, v'^^ = 2.552, v'^^ = 0.0074, 
and f^'g = 11.83 we obtain the masses: m^ = 0.006, mg = 0.1, = 4.2. 
For the w-type quarks we use, a„ = a^, but now hu = h'^ = 3, and v[ = 
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53.8, v'^ 



0.4, 



0.0001 and v'^_ 



0.1, v',^ 



0.0025, mc = 1.27, and mt 



124, v',^ 



0.082, which give 
171.2. With these VEVs and parame- 
ters we obtain V^'^ matrices given, according with the definition above, the 
mixing matrix 



^ 0.97747 0.21105 0.00148^ 



\V( 



CKM\ 



0.21083 0.97673 0.03923 
0.00972 0.03803 0.99923 



'111 



Comparing the elements of the matrix above with experimental data in Sec- 
tion 11.2 of PDG [2S], we see that all the matrix elements, but us and uh 
entries, are within the level of confidence (1-2) a. However, this is enough to 
show that a realistic mixing matrix may arise from the quark mass matrices 
given in Eqs. ([6]) and ([7j). A more rigorous and detailed analysis is being done 
and will be presented later. 

The suppression of the FCNCs coupled to neutral Higgs bosons can be ob- 
tained not naturally in the sense of Ref. [26], but at least as a reasonable fine 
tuning among the mixing matrix in the neutral scalar sector, f/j^, or by the 
mass of heavy Higgs scalars. For instance, consider the neutral kaons mass dif- 
ference, ^Mk ~ C./i^i.(E, mV^]) where = K I )adP , 

with symmetry eigenstates, say iff {i fixed), and physical scalars hj with mass 
TTij. In order to be consistent with data C,l^ must be smaller than the contri- 
bution of the SM: i.e., smaller that Glml{V*aV,,f / IQii"^ = 0.5 x lO-^GeV-^ 
(for simplicity we have taken into account only the c quark contribution and 
neglect the QCD and m^/M^ corrections). Notice that the sum over the con- 
tributions of each i is made coherently. The expression above is valid for any 
scalar heavier than niK- For instance, h'^\J2a(yL)as(yR)ad\'^ ~ 10~^ ~ 10~^ 
and |?7jjp/m| ~ (10^^ — 10^^)GeV^^ we have the appropriate FCNC suppres- 
sion. These constraints are not a large fine tunings since it can be satisfied 
with a large set of unitary matrices and/or the neutral scalar mass spectrum 
and mixing. Recall that there is no more freedom to redefine the quark phases 
since they all have already been used to put the CKM matrix in the usual 
form, hence has five physical phases. The case of FCNCs in the u-quark 
sector involve the matrices and Vj^ . 

On the other hand, it would be interesting if we can show that appropriate 
hierarchies among the VEVs arises from an analysis of the scalar potential j27j. 
We have verified that, under reasonable assumptions on the parameters in the 
scalar potential, this is indeed the case. For the sake of simplicity we will 
consider the scalar potential involving only the scalar doublets that couple to 
quarks: the triplets H', $' for w-like quarks, and H", $" for d-\ike quarks. 
Denoting X„ = H'^H' and = $'t$'^ = H"^H" and Yd = we 
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write the most general scalar potential invariant under the SM and ® ^3 ® 
Z3 Z3 symmetries as: 

+ 9'AH"emH"emt + ^{{H' e^"UH" e^%) + i/.c). (12) 



We have omitted summation symbols and used the notation p, p' = u, d; mn = 
11, 1'l", AA; kl = 11, 1'l", BB; a/3 = 11, 1'l", AB; and st = 11, 1'l", AA, BB, AB. 
For example, XPPi[X,UX,,l = A?f[Xp]i[Xp,]i+A?f;,[Xp]i4Xp,]i»+A^^;[Xp]^[Xp,]A 
+ ■ ■ ■, and each term having contributions p, p' = u, d. Considering the Eq. f|T2l) 
we obtain the constraint equations arising from the conditions dV/dv'/ = 0, . . .. 
For simplicity, in the following analysis it will be also assumed that the VEVs 
are all real. In this case the condition 1 + u + uj'^ = implies that in the 
constraint equations the coupling constants will appear as the sum of some of 
the coupling constants in Eq. (fT2|) . The inclusion of all the scalar fields of the 
model and also the constraints coming from the condition d'^V / dv'l dv'- > 0, . . . 
will be shown elsewhere. We have obtained the constraint equations under the 
following conditions: v'l = v'^^, Vz; and A^„ = K^f, Wmn; kl; and verified that, 
among the solutions to those equations, it is possible to have the following 
hierarchy: ^ fg ^ vi^. However for generating the quark masses, as we 
have shown above, it is necessary that v'}"'^ ^ '^'j"'^ ■ Notwithstanding, it is 
interesting that such hierarchies can arise from the model, and we hope that 
a more detailed treatment of the scalar potential will produce a more realistic 
hierarchies among the VEVs. 

As in any model with broken discrete symmetries there are potential troubles 
with domain walls [2H]- However, this issue can be solved by allowing soft 
terms in the scalar potential violating these symmetries. 

In conclusion, the mass matrices obtained, Eqs. ([5]), ([6]), ([7]), and ([8]), arising 
because of the symmetries of the model, give appropriate insight concerning 
the solution of the flavor problem generating at the same time a realistic 
VcKM matrix. Of course, it is necessary to explain how these symmetries are 
realized from a more fundamental theory. Our scheme has some similarities 
with the "private Higgs" of Porto and Zee [29] in which there is one Higgs per 
fermion and a general mass matrices. However, we have more than one Higgs 
per fermion and the mass matrices have a specific texture of zeros. 

We thank an anonymous referee for useful comments. ACBM was fully sup- 
ported and VP partially supported by CNPq. 
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